BROWN-HET-996 
hep-th/9506063 

Quenched Two Dimensional Supersymmetric 

Yang-Mills Theory 

Miao Li 

Department of Physics 

Brown University 
Providence, RI 02912 
li@het . brown . edu 

By studying the pure Yang-Mills theory on a circle, as well as an adjoint scalar coupled 
to the gauge field on a circle, we propose a quenching prescription in which the combination 
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1. Introduction 



The two dimensional supersymmetric Yang-Mills theory may turn out to be a two di- 
mensional matrix model which can be solved in the large N limit. A set of loop equations 
are solved in with a certain assumption, its spectrum is studied numerically in in 
the light-cone formalism. There seems to be at least two major motivations for studying 
supersymmetric gauge theories in various dimensions. The first is the hope that deeper 
understanding of these theories may shed light on several longstanding open problems in 
particle theory and quantum field theories, such as confinement [Q, dynamic supersymme- 
try breaking [|J and various kinds of duality P|| . The second motivation, which is equally 
important if not more important, is the search for higher dimensional solvable matrix mod- 
els. As we have learned from the study of zero and one dimensional matrix models, string 
theory can be reformulated in terms of matrix models. In such models, some spacetime 
dimensions are dynamically generated, and powerful mathematical tools can be developed 
to calculate various physical quantities. It may be justified to hope that string theory will 
be ultimately formulated in this fashion. This latter motivation serves as the primary one 
for our previous study £j] and the work presented in this paper. Recently, based on the 
explicit solution of the low energy effective action in the supersymmetric SU(N) Yang- 
Mills theory in four dimensions |J, Douglas and Shenker extract some interesting large N 
information about the theory [[7], and find that in the large N limit, the validity region 
of the low energy solution becomes very narrow. Because of the importance of the large 
N solution of super-gauge theories, it is then desirable to explore all possible valuable 
methods. 

As long as the large N problem is concerned, Eguchi and Kawai showed that the 
Euclidean lattice gauge theory can be reduced to a model at a single site with D matrices 
H , here D is the spacetime dimension. Unfortunately, it was shown subsequently that the 
EK model suffers breakdown of global U(l) symmetries at weak coupling in dimensions 
higher than two 0, so the EK model is not capable of reproducing the large N result 
of the Wilson theory in the physical weak coupling regime. A quenched model was then 
introduced in ||. This idea was substantially augmented by Parisi and by Gross and 
Kitazawa JllJ, these authors were able to show that the quenched matrix model indeed 
correctly produces the large N results for correlation functions in the weak coupling regime. 



We shall follow refs.[]T2| to work with the Hamiltonian formalism. 

The quenching prescription proposed in this paper is slightly different from the conven- 
tional one. Instead of first calculating quantities with fixed quenching momentum matrix, 
and then averaging over momenta, we propose to consider the combination of the mo- 
mentum matrix with the spatial component of the gauge field as a dynamic variable. The 
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motivation for doing this is from consideration of the pure Yang-Mills theory on a cylinder, 
in which the theory already reduces to a quantum mechanics problem. By construction 
this prescription works only in two dimensions. The quenched super gauge theory, being 
equivalent to the original theory in the large N limit, is an one dimensional N = 1 super- 
symmetric multi-matrix model. There are two Hermitian bosonic matrices, two Hermitian 
fermionic matrices, or equivalently one complex bosonic matrix and one complex fermionic 
matrix. This matrix model is interesting in its own right. It is interesting to study its 1/N 
corrections, even these may have nothing to do with 1/N corrections in the original super 
gauge theory. In the similar spirit, one may try to generalize this super symmetric matrix 
model in various ways, and to ask the question whether a doubling scaling limit exists, 
and if it does, what kind of string theory it describes. 

We shall not try to solve the quenched model directly in this paper. The main result of 
this paper perhaps is the construction of a solvable model by extending the supersymmetric 
quenched model to include one more complex fermionic matrix. There is still one complex 
bosonic matrix. The bosonic part of the action differs from the N = 1 model by a potential 
term. The N = 1 supersymmetry in the extended model can be explained as Parisi-Sourlas 



supersymmetry [[13|. Since such a model can be reduced to a Gaussian model together 
with a pair of stochastic equations, Green's functions of bosonic matrices are calculable in 
principle, once the first order stochastic equations are solved. It is conceivable that further 
understanding of the extended model will shed light on understanding of the quenched 
model. 

Since the quenching prescription was proposed more than ten years ago, it is appropri- 
ate to review the essential ingredients first. We shall do this in the next section. We then 
proceed in sect. 3 to discuss the prescription for two dimensional gauge theories by starting 
with the pure Yang-Mills theory defined on a cylinder. We shall argue that by promoting 
the combination of momentum matrix and the spatial component of the gauge field to a 
dynamic matrix, the average over momenta is automatically done. Also, it is necessary to 
rescale the coupling constant. The consistency of this prescription is checked in the case 
of an adjoint matter coupled to the gauge field. The supersymmetric quenched model is 
introduced in sect. 4. We do not know how to solve this model yet. Then in sect. 5 we 
extend this model to obtain an exactly solvable model. This model exhibits Parisi-Sourlas 
supersymmetry, and associated first order stochastic equations can be integrated. The last 
section is devoted to a discussion. 
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2. A Brief Review of the Quenching Prescription 



To illustrate the idea of quenching, consider a two dimensional Hermitian matrix 
M(x,t) with the following action 

S(M) = J d 2 xtr Q(d t M) 2 - \{d x M) 2 + -^=M s ^j , (2.1) 

where the coupling constant g is held fixed in the limit iV — > oo. 

Instead of quenching all spacetime momenta as in ||, 10] , lj| , only the spatial mo- 
mentum is to be quenched in this paper [TJ]. The prescription is to replace the spatial 
dependent matrix field M(x : t) by exp(iPx)M(t) exp(-iPx), here P is a diagonal matrix 
with eigenvalues p^. The derivative d x M is replaced by i[P, M] and the quenched action 
reads 

S = a J dttr Q(<9 t M) 2 + ~[P, M] 2 + ~^ m3 ) > ( 2 - 2 ) 



where a = 2n / A is the ultraviolet cut-off whose utility we will see shortly, M in the above 
action depends only on t. The propagator derived from (|2.2|) is 



(M ij (t)Mi k (t')) = -~8 ik 8 jl (d?+p 2 j )- 1 8(t-t'), (2.3) 



a - lJ/ 

where pij = Pi — Pj, the momentum carried by Afy. We thus see that this propagator is 
almost the same as (Mij(p, t)Mik(q, t')) in the original model if one lets p = Pij, except for 
lacking of the delta function factor S(p + q). 

We now show that any planar vacuum diagram in the original model is recovered by 
averaging the corresponding diagram in the quenched model with an integral 

•-fir* ^ 

fJi J -A/2 2tt 

each integral in the above product is normalized to 1. To see this, consider a planar diagram 
with I loops, n vertices and p propagators. With the standard double-line representation, 
on assigns momentum pi the line with index i, then momentum conservation is automatic 
in the quenched model. An example of vacuum diagrams is drawn in the following figure. 




The propagator is essentially the same as in the original model, with an extra factor 
a -1 , so there is an additional factor a~ p from all propagators. The n vertices contribute 
a factor (ag) n N~ n / 2 and the contraction of matrix indices gives rise to a factor N l+1 . 
So altogether there is a factor a n ~ p g n N l+1 ~ n / 2 = a n ~ p g n N 2 , where relations 3n = 2p 
and / — p + n = 1 are used. The factor iV 2 is the right one for a planar diagram. Now 
multiplying the result by the integral factor ( |2.4[ ), I integrals in this factor are identified 
with loop integrals in the original diagram, leaving a factor a 1 together a -1 in ( |2.4| ). The 
rest integral factors normalize to 1. Thus, the a dependent factor is a l ~ p+n ~ 1 = 1, with 
the planar relation I — p + n = 1. This explains why a factor a in the quenched action is 
introduced. In the limit N — > oo all planar diagrams are included, thus the free energy 
calculated with the quenching prescription is the same as in the original model. 

Calculating Green's functions requires a little modification. First, according to the 
quenching prescription, tr M n (x,t) = tr M n (t), so Green's functions of these quantities 
will be independent of positions. This is true in the large iV limit according to the factor- 
ization theorem. Next, one would like to calculate the expectation value of quantities such 
as tr [M(x, t)M(y, t')] or tr [M(p, t)M(q, £')]. By the quenching prescription 

(tr [M(p,t)M(q,t')}) = (2tt) 2 5( P - Pij )8(q + Pij ) (M(t) ii M(t / )^), 

where sum over indices is assumed. In the leading order, applying ( |2.3|) and ( |2.4|) , one 
finds that the result is the same as in the original model except for an additional factor 
a -1 . It is easy to see that this factor exists for all connected planar diagrams contributing 
to this Green's function. Similarly, when one considers a connected Green's function of 
this type involving n matrices, an additional factor a -1 is to be compensated. 
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3. Quenched Yang-Mills Theory on a Cylinder 

3.1. Pure Yang-Mills Theory on a Cylinder 

The pure Yang-Mills theory on a cylinder reduces to a quantum mechanics problem 
I(|[l7|], since the only dynamical degrees of freedom are winding modes of the gauge field. 



Because the system itself already effectively collapses to a "point" before quenching, it is 
desirable to compare what obtained in the quenched model to the un-quenched model. 

We start by showing how the pure Yang-Mills theory effectively reduces to a quantum 
mechanical system. The standard action 

S = -^ [ d 2 xtr (F 01 ) 2 (3.1) 



2g 2 j 

can be cast into a first order form, upon introducing the canonical momentum II D 

s = Jd 2 xtr (-^(n D ) 2 + n^F 01 ) 

= J d 2 xtr (-^(n°) 2 + U D d t A + A D x U D ^j , 



(3.2) 



where A = A\. Integrating out A in the path integral f [dA dAdU D ] exp(iS) results in a 
delta function 5(D X U D ). This forces U D to satisfy D X U D = 0. The solution is given by 

n D (x) = U- 1 {x)U. D U{x\ U(x) = Pexp(i [ Adx), 

Jo 

where U D is n D (0), a matrix independent of x. Let the circumference of the circle be L. 
By the periodic condition, U D (L) = II D (0) = U D , one finds [U, U D ] = 0. Here U is the 
holonomy around the circle U = U (L) . Thus, upon substituting this solution into the path 
integral, the delta function reduces to 6([U, H D ]) and the path integral itself becomes 

[dAdn D ]5([U,n D ])exp(i j dttr (~^(T1 D ) 2 + T1 D J dxU(x)d t A(x)U- 1 (x))^j . 

It is easy to see that J dxll (x)dtA(x)U~ 1 (x)) = —idtUU , from the definition of U. 
Thus, all degrees of freedom in A(x) collapse to that of U, and the above path integral 
finally reduces to a path integral of a quantum mechanic system 



J [dUdU D ]5([U,U D ])exp (i j dttr {-^{H D f - m D d t UU~ x ] 



(3.3) 



where the appropriate measure [dU] is the Haar measure. Instead of working with the 
holonomy, one can work with D with the definition U = exp(iDL). Using the following 
formula 

-idtUU' 1 = L I dT( lL ' LT U t De- 
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in ( |3.3|) and noting that exp(iDLr) effectively commutes with Il D , thanks to the delta 
function in the path integral, one finds that these exponentials cancel upon taking trace. 
The delta function 8([U, H D }) can be replaced by another one S([D, H D }) and finally 

J [dDdU D ]S([D, T1 D ]) exp (i J dtLtv (-^-(U D ) 2 + U D d t D)^j , (3.4) 

where extra care need be exercised in definition of the measure [dD]. One can define the 
measure as the usual flat one, and the nontrivial factor coming from the Haar measure 
[dU] can be absorbed into the definition of the delta function. Eq. (|3.4| ) is our final formula 
for comparison with the quenched model. 

Following the prescription given in the previous section, we use the substitution 

U D (x) = e iPx U D e - iPx , Aq{x) = e lPx 'A e~ lPx , A(x) = e lPx Ae~ lPx (3.5) 

with a diagonal matrix P whose entries take the form pi = 2imi/L, m is an integer. 
Furthermore, D X TL D = i[D,H D ], where D = P + A. We shall see presently that D is to 
be identified with the previously introduced D in ( |3.4[ ). The action ( |3.2| ) is replaced by 

S = aJ dttv ^-^(n D ) 2 + U D d t A + iA [D, U D ]^ . (3.6) 

Again integration of Aq results in constraint [D, U D ] = 0, and the path integral becomes 

J [dDdU D ]S([D, U D ] exp (i J dtatv (-^-(U D ) 2 + U D d t A)^j . (3.7) 

To see what quantity corresponds to D, consider the Wilson line U xy = Pexp JJ A(x)dx. 
With the quenching substitution ( |3.5|) , U xy = exp(iPx) exp(iD(y — x))exp(iPy), and in 
particular for the Wilson line which wraps the circle once U = exp(iDL). We thus see that 
D is the same quantity as introduced in ( |3.4| ). With this in mind, the path integral ( p.7|) is 
seen to be almost the same as the un-quenched version ( |3.4[ ). The differences between the 
two are that it is A appears in the action of the quenched version, and that the ultraviolet 
cut-off a appears in the quenched action while it is the infrared cut-off L appears in the 
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un-quenched action. The first difference is easily resolved by replacing dtA by dfD, since 
P is supposed to be held fixed with discrete eigenvalues. 

Before resolving the second difference, we remind ourselves that in the original pro- 



posal of 1 11], the constraint that the eigenvalues of D coincide with those of P is imposed. 



This amounts to inserting the following factor 

J [dU}5(D - UPU-^Aipi), A(pi) = HiPi-Pj) 2 (3.8) 

i<j 

into the path integral. We now argue that such constraint should not be imposed in 
our case of a compact circle. As we have seen, an eigenvalue of DL is the phase of an 
eigenvalue of the holonomy U. With the above constraint, this eigenvalue would be 27raj 
and the corresponding phase is trivial in the holonomy. So the constraint ( |3.8|) would lead 
to a trivial theory. 

Now the second difference mentioned before between the two actions is formally re- 
solved by rescaling U D — > j^n D and g 2 — > j^g 2 - The rescaling of H D does not change 
the quenched model, since it is just a matter of convention. The rescaling of g 2 is a little 
disturbing. By viewing D, instead of A, as the dynamic degrees of freedom, this rescaling 
appears necessary in order to recover the original theory at the large iV limit. Such mod- 
ification is not absurd as it might appear. Note that except for the Gauss law constraint 
[D, U D ] = 0, the action in (|3.7[) is quadratic, so a straightforward perturbative argument 
as presented in the previous section is lacking. In particular, averaging over pi with ( |2.4| ) 
can not be introduced directly. Finally, we would like to point out that the ratio L/a can 
be taken equal to N. The reason is the following. The number of possible values of pi with 
both a ultraviolet cut-off and an infrared cut-off is just L/a. When P is absorbed into D, 
the rank of D, N, counts effectively the number of possible values of pi. This is pointed 



out in [[14] as well as in ||15|| . 

To check the consistency of our special prescription, we shall consider the system of 
Yang-Mills coupled to an adjoint scalar field in the next subsection. 

3.2. Yang-Mills Coupled to an Adjoint Scalar 

Let <p be an adjoint scalar, therefore a Hermitian matrix field for a gauge group U (N). 
In this subsection we are interested in studying the coupled system with an action 

S = ^2 / d2xtr (( F oi) 2 " (D^) 2 ) , (3.9) 



with definition of covariant derivatives D^cj) = d^4)+i[A^, <f>]. In the above we appropriately 
rescaled <p so that the action is weighted by a factor 1/g 2 . This action can also be written 
in a first order form 

s = j d 2 xti (-^[(n D ) 2 + (n*f} + n D F 01 + n*D t( i>-^(D x( t>)A 

= Jd 2 xtr ^-^[(U D ) 2 + (U^) 2 ] + U D d t A + U^d t( p-^(D x( p) 2 ^ (3.10) 

+ J d 2 xtr A (D X U D + 

Again integrating out Aq in the path integral imposes the Gauss law D X U D + i[(f>, EP] = 0, 
which is solved by the following expression 

U D (x) = U- 1 (x)U D U(x)-i [ dyU(x,y)[4>(y),U%)]U(y,x), (3.11) 

Jo 

where H D = U D (0), and U(x) is defined as in the last subsection, and U(x,y) = 
_Pexp(i A(x)dx). Now if one puts the system on a circle of circumference L, the periodic 
boundary condition for TL D (x) no longer results in [U, Il D ] = 0, because of the second term 
in ( |3.11|) . Also, after substituting ( |3.11|) into the action (|3.10|) , the action not only depends 
on U = U (0, L), but also on U (x, y). Therefore, in the quenched model, the eigenvalues of 
.D, as defined in the previous subsection is no longer restricted to a circle, since U (x, y) is 
not well defined on this circle. Thus, unlike in the pure Yang-Mills theory, the momentum 
as part of an eigenvalue of D will enter in the story in an essential way. 

The solution ( |3.11| ) of the constraint can be simplified by using gauge transformation 

(f>(x) -> U-\x)(j)(x)U(x), n^(x)) -> U- 1 {x)Il^{x))U{x) 1 

Tl D (x) -> U- l {x)Il D {x)U[x), A(x) -> U- 1 (x)A(x)U(x) - iU- 1 (x)d x U(x). 



The action (|3.10|) is invariant under this transformation, since the additional term resulting 
from the transformation is proportional to the constraint, therefore vanishes. The solution 
( [3. 11|) after subject to this transformation reads 

Tl D (x)=n D -i f X dy[<l>(y),IL*(y)]. (3.12) 

Since U (x) is not necessarily periodic, the periodic boundary conditions for new fields are 
twisted by U, for example H D (L) = L r Il r) (0)L r_1 . With the help of the above solution, 
this condition translates into 

[U,Tl D } = -i [ L dy[<P(y),U^y)]]. (3.13) 
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One can substitute ( |3.12|) into action ( 3.10 ). Once again, one would obtain a term 

9 2 L 



J dttv (n D ) 2 



in which the infrared cut-off L appears explicitly. This suggests that the quenching pre- 
scription introduced in the previous subsection also works in this model. To further confirm 
this, we suggest an unconventional perturbation scheme in which the whole term 



1 



2r 



d 2 xtr (D a 



is treated as a perturbation term. Certainly there are other terms resulting from solving 
the constraint. Let us focus on the above term. To normalize the kinetic terms in ( |3.10|) , 
one has to rescale (f) — > gcf>, — > TL^/g, D x — > gD x . Thus, the above perturbation becomes 

-y / <4l>^ 2 ' (3 - 14) 

n 

where we replaced the integral over the spatial dimension by a sum over discrete momentum 
modes. Accordingly, the quantity tr (D x (p) 2 should be viewed as the Fourier transform of 
the original term. Essentially, there is no ultraviolet cut-off involved in this perturbation 
theory. D x (f) is effectively replaced by iptfi + i[a,(f)] where a is the constant mode of the 
gauge field. For a detailed discussion we refer to [[19] where an adjoint fermion coupled to 
the gauge field is discussed. Now one can replace the sum over momentum modes in ( |3.14| ) 
by a trace, thus the perturbation ( |3.14| ) becomes a four point vertex and the combination 
of p + A becomes an independent field. 

Finally, we are in a position to propose our quenching prescription. First, as usual, 



A(x) = e lPx Ae~ lPx , U D (x) = e lPx Ii D e~ lPx , 
(j)( x ) = e lPx (pe- lPx : n*(x) = e lPx W t> (x)e- iP 

With these substitutes, the action in ( |3.10|) becomes 



(3.15) 



S = L J dttr (- 9 l[(Jl D f + (U^) 2 }+U D d t D + U^d t( p+ J-^D,^, (3.16) 



with constraint 



[d, n D ] + [cf), n^] = o. (3.17) 



It is easy to see that the above constraint is a consequence of the old one with prescription 
( |3.15|) . What is new here is that the usual ultraviolet cut-off a is replaced by the infrared 
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cut-off L, and D is promoted to a dynamic variable so the average over momentum is not 
necessary. This is also why a ultraviolet cut-off is not necessary, since the propagator of 
4> does not involve P. The necessity of an infrared cut-off is best shown in study of the 
pure Yang-Mills theory To see that this prescription is consistent, we note that any gauge 
invariant quantity, such as 

tr [<p(x)U (x, y)<fi(y)U(y, x)] = tr 4>exp(iD(y — x))4>exp(iD(x — y)), 

contains P through D. The prescription works only for gauge theories, since in a matrix 
model with only global symmetry, P appears explicit in all correlation functions. 

Before closing this section, we mention once again that the Gross-Kitazawa constraint 
(|3.8|) should not be imposed in this model on a compact circle too. And the eigenvalues 
of D, unlike in the pure gauge theory, should not be restricted to live on a circle. 

4. The Quenched Supersymmetric Yang-Mills Theory 

4-1. The Quenched Model 

In this section, we shall introduce the quenched supersymmetric gauge theory. First, 
let us introduce the 2D super-gauge theory. The vector super-multiplet in two dimensions 
consists of a vector field, a scalar field and a Majorana fermion. The scalar field is needed, 
since the vector field in two dimensions does not have dynamical degrees of freedom. The 
action 

S=^J d 2 xtv (-^-^{D^-iX^D^X-Xa^X]^, (4.1) 

is invariant under the SUSY transformation: 

SAfj, = -2zAa M e, 
8(j) = 2i\a 3 e, (4.2) 
5X = o^Foie - a^a 3 D^e. 

We use notations of [p!8| , o" M = (1,— a 1 ). The constant g 2 N is held fixed in the limit 
iV — > 00. Here for simplicity we consider the gauge group U(N). 

It is necessary to fix a gauge in the Hamiltonian formalism. An arbitrary gauge will 
spoil supersymmetry, thus many advantages brought about by SUSY may disappear. As 
was shown in there is a gauge in which half of SUSY is broken, but half of SUSY 
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survives. In this gauge, Aq + <p = 0. Thus we are left with two bosonic fields A = A\ and 
</>, and two fermionic fields A a . The conjugate momenta are 



U A = \F i = \(d A + D x c 

g 2 g 2 

= -1^00, n A = 
g g 



(4.3) 



with the standard commutation relations 

\^-ij{ x )-> A ik{y)] = -i8i k 5ji8(x-y), [ut j (x) 1 (pik(y)] = -i5 ik 5ji5(x - y), 

••i* g 2 (4 - 4) 

{X^{x)^{y)} = g -8 ik 8 jl 8 aP 8{x-y). 
Now the Hamiltonian and the unbroken super-charge Q are, respectively 

H = g 2 J dxtr QtIT 4 ) 2 + ^(IT^ - J dxtr U A D x( f> 
+ -1 J dxtr Q(^) 2 -zAa 1 D x A + A((T 3 -l)[(/),A] 

Q = J dxtr ^AiII* + A 2 (IT A - ^D x </>fj . 



(4.5) 



and 



(4.6) 



The relation between these two quantities is H — {Q, Q}. 

Now we are in a position to introduce the quenched model. The model is obtained 
by replacing every matrix field F(x) with exp(iPx)F exp(-iPx), where F is independent 
of x and P is a real diagonal matrix, as in the previous section. The covariant derivative 
D x F(x) is then replaced by exp(iPx)i[D, F] exp(— iPx), where D = P + A. On a circle, 
the eigen- values of P take discrete values 2-rm/L, n is an integer. With these substitutions, 
the Hamiltonian and the super-charge are truncated to a point 

H=^tr ((II D ) 2 + (II*) 2 ) -zLtr II D 

L ( 1 \ (47) 

+ -tr ( — -[-D, 4>] 2 + Xa l [D, A] + A(a 3 - 1)[0,A]J . 



and 



Q = Ltr ( X 1 U* + X 2 (U D - ^i[D, </>]) ) , (4.* 



where we replaced II by II D . Still these quantities are not convenient to work with, since 
the commutation relations are given by (fOJ) with 5(x — y) replaced by 1/L. Note that here 
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(4.11) 



our prescription is different from the usuai one, in that the delta function is not regularized 
by 1/a with a short-distance cut-off a. The argument is the same as we presented in the 
last section. The commutators are simplified by substitutions 

D ^7z D > ^7z^ x ^7z K 

u D -> —n D -> — ^ir^ 

Explicitly 

[n.?,Dik] = -i5i k 5ji, [U.fp(j)ik] = —iSikSji, 

Ik 1 ^ 

The super-charge reads 

Q = tr (Ain^ + A 2 (n D -zZ[D, 0])), (4.10) 
where the parameter / = g/y/L. The final ingredient is the gauge transformation 

D^UDU' 1 , (p^U^U' 1 , X^UXU' 1 , 

where the original gauge transformation parameter U(x) is also replaced by 

exp(iPx)U exp(-iPx). 

The above gauge transformation is generated by the generator 

G = i[D,U D ] + 2A 2 . (4.12) 

A physical state is annihilated by G. 

4-2. The Superspace Formulation 

The super-charge ( |4.10|) can be written in a more symmetric fashion, if one introduces 

A = Ai — zA2, A = Ai + zA2, 
M = (p-iD, M = <p + iD, ( 4 _ 13 ) 

A = n M — ^ [M, M] , I = ^[M, Ml 

then 

Q = tr (AA + XA), (4.14) 
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with corresponding action 

S = J dttr [^d t Md t M - j[M, M}d t (M - M) + iXd t X 
- l -\(M, X] - l -X[M, X] + l -X[M + M, A]]. 
This action is invariant under the following SUSY transformation 

5M = -2ieX, 6M = -2ie\ 
SX = ed t M, 5X = ed t M. 



(4.15) 



(4.16) 



Since the above transformation has a simple form, one would guess there might be a 
simple superspace formulation. Indeed there is. Introduce a real fermionic coordinate 9 
as the superpartner of t. Let D = dg — i6d t be the super-covariant derivative. The super- 
charge commuting with D is Q = dg + i6dt. There is {Q, Q} = 2idt = 2H. Introduce the 
following super-fields 

$ = M + ^A, $ = M + i6X. (4.17) 



It is easy to see that the free part of the action (|4.15|) is 



-i J dtdOtr D$D 2 § = J dttr i (d t Md t M + iXd t X) , 

where there is an extra factor 1/2 for the fermionic part, which is due to a rescaling 
A — > X/y/2. The cubic terms in ( [4 . 1 5| ) can also be sorted out easily. We simply write down 
the whole action in terms of the super-fields 

S = J dtdOtr ^-^D$D 2 $+ ^ . (4.18) 



5. An Extended Model With Parisi-Sourlas Supersymmetry 

We have tried to solve the quenched model introduced in the last section, for example, 
by constructing some conserved quantities. So far we have not succeeded in tackling this 
model directly. The reason why any simple minded method of constructing conserved 
quantities does not work is that all A and A introduced in ( |4.13|) do not commute. For 
example, a commutator of two matrix elements of A depends on M. 

In this section, we suggest to study an exactly solvable model with one more complex 
fermion. This model is very similar to the N = 1 model introduced in the previous section. 
We believe that better understanding of this model shall shed light to the N = 1 quenched 
model. 
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The supersymmetry to be introduced here is Parisi-Sourlas supersymmetry. Therefore 



it is convenient to follow the line of the original papers |T3[ to introduce this model. We 
start with an action obtained by adding a potential term 

^tr [M,M] 2 , 



to the bosonic part of ( |4. 15 ) . This term is positive definite, since \M,M] is Hermitian. 



This means that a negative definite term is added to the quenched Hamiltonian. The so 
obtained bosonic action is 

S(M,M) = j dttv[^(d t M-^[M,M])(d t M + ^[M,M])]. (5.1) 

Note that this action is a complete square. 
Define the complex Gaussian fields 

v = d t M -^[M,M], f ] = d t M + ^[M,M], (5.2) 



JidndfjleU^^. (5.3) 



with the measure 



The Green's functions of M and M are defined through solution to ( |5.2| ) together with the 
above Gaussian path integral. Denote the solution (with appropriate initial conditions) of 
( |01) by M = M(tj, fj),M = M(t], fj). A Green's function is then 



[dr]df})F[M{r], fj),M(r],fj}eif dttT m 
= J[dMdM}det[^^]F[M,My s : 

where in the second line we changed the integral variables from the 77's to the M's, and S 



is the action given by (|5.1|). The Jacobian, following relations (|5.2|), is given by 

(5.4) 



d(V:V) _f9t-f[M,.] -f[-,M] 



d(M,M) V -fh^] dt~%[M, 
Its determinant can be expressed as a fermionic path integral. Let ^ = [A, ijj] , then 
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with 

W T) = */*t, W ,#)^(A,X)< 

= J dttr (i^d t X + iifjd t X + ^(i(j[M, X] + $[M, A] - 4>[M, A] - $[M, A])^ . 

(5.5) 

Note that the structure of this action is almost the same as the fermionic part of action 
( |4.15|) , except here there are two complex fermionic matrices. (Formally if ip — > X and 
ij) —>■ A, then the above action is the same as in ( f4.15|) .) Now a Green's function of M and 
M can be written as 

(F[M,M}) = J [dMMd^d^]F[M,M]e l[s{M;M)+s ^^ )] . (5.6) 

As always, the combined action S(M, M) + S(^, ^) possesses supersymmetry, Parisi- 

Sourlas supersymmetry. Unlike in the simplest cases, the supersymmetry parameter is a 

real Grassmanian: 

SM = -2ieX, SM = -2ieA, 

difj = e(d t M + ^[M, MM , 



( 5 - 7 ) 

8i; = e(d t M- i j[M, M] 

X and A are invariant under this supersymmetry. We suspect that there is one more 
supersymmetry, perhaps a nonlinear one, but have not been able to sort it out. The 
supersymmetry ( |5.7|) itself does not generate the Hamiltonian, precisely because the A's 
are invariant under it. To be able to generate the Hamiltonian, one perhaps need another 
supersymmetry. 

The first order nonlinear differential equations in (|53) are integrable. Observe that 
the "friction" term il/2[M, M] is anti-Hermitian, so the Hermitian part of M satisfies a 
linear equation. Use the same notation as in (|4.13|), denote this Hermitian part by <p and 



the anti-Hermitian part by D. A general solution for </> is 

(p(t)=f dt'Re-n(t')+(/)(-oo), (5.8) 

J —oo 

where Rery denotes the Hermitian part of 77, 0(— 00) is the value of (p at the infinite past. 
Now D satisfies the following equation 

d t D -il[(j),D] = -Imrj, (5.9) 



15 



with Imi] being the anti-Hermitian part of r\. Since <p(t) is known as in ( |5.8| ), the solution 
to the above equation can be expressed in terms of </>, 

D(t) = - f dt'U- x {t',t)lm.ri{t')U{t',t) +D(-oo), (5.10) 



where 



U(t',t)=Pexp(-il [ <p(r)dr 
Jt' 



That (p appears in the above path ordered integral is not surprising. Recall that in the 
quenched model, <p is J us t — ^-o in the gauge in which one supersymmetry is preserved. 
Thus, U(t' , t) is just the Wilson line along the time direction. 

Now any Green's function involving only matrices <p an d D, or equivalently M and 
M, can be calculated using (|5.8| ) and ( |5.10| ) together with the Gaussian path integral (|5.3|). 
We shall not try to develop a systematic technique to do this in the present paper, but 
shall leave it for future work. Since the 77's are Gaussian fields, it is possible to treat the 
large iV problem with the help of the master field, for these are the so-called free variables. 
Some Green's functions involving fermionic matrices can also be calculated, using the Ward 
identities associated with SUSY in (]5.7|). For instance, starting with 



(8 e tv [V>M]) = 

one derives 

(tr [tp\]) = l -(ti [M(d t M + ^{M,M])}). (5.11) 

This Ward identity holds for finite N. There is no anomalous term, since SUSY is not 
spontaneously broken. This is guaranteed by the fact that the solution to the Langevin 
equations ( |5.2[ ) is unique for given initial data. 

Finally, a technical issue to be addressed is the projection to the singlet sector. The 
solution is not difficult and again we leave it to future work. 



6. Discussion 

We have shown that the two dimensional supersymmetric Yang-Mills theory truncates 
to an one dimensional supersymmetric matrix model in the large N limit. Although the 
quenched model may have nothing to do with the original SYM theory beyond the large 
N limit, it is tempting to treat itself as an interesting matrix model, and to ask various 
questions as often arise in dealing with a matrix model. Before the large iV problem 
becomes tractable, asking these questions remain un-practical. 
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The extended matrix model studied in sect. 5, exhibiting Parisi-Sourlas supersymme- 
try, can be solved exactly. This model is very similar to the quenched model, so further 
study will undoubtedly shed light on understanding of the quenched model. Some mech- 
anism is needed to truncate the number of fermionic matrices, in order to come down to 
the quenched model. We have studied the quenched N = 2 two dimensional SYM theory, 
and the resulting theory contains the same number of fermions as in the extended model. 
There is N = 2 supersymmetry and an additional bosonic Hermitian matrix. We intend 
to publish some results concerning this model elsewhere. For now, we just mention that 
the N — 2 2D SYM theory is the dimensional reduction of the N = 1 4D SYM theory, and 
the latter is shown to possess the Nicolai mapping [ 20j . In a sense, the Nicolai mapping is 
nothing but a generalized Parisi-Sourlas mapping. Therefore, we suspect that this model 
might be interconnected with the two models studied in this paper. 

It is our hope that some higher dimensional string theory will eventually emerge 
from these supersymmetric matrix models. With multi-matrices at hand, it is possible 
to generate more than one dimensions. Also, so far only supersymmetry can help in 
overcoming the so-called c = 1 barrier, and in eliminating problems such as tachyon. We 
hope to return to models introduced in this paper in near future. 
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